REFLECTIONS ON RAMANUJAN'S 
MATHEMATICAL GEMS 



G. D. ANDERSON AND M. VUORINEN 

Abstract. The authors provide a survey of certain aspects of their 
joint work with the late M. K. Vamanamurthy. Most of the results are 
simple to state and deal with special functions, a topic of research where 
S. Ramanujan's contributions are well-known landmarks. The compre- 
hensive bibliography includes references to the latest contributions to 
this field. 



1. Introduction 

1.1. Ramanujan's life. Srinivasa Ramanujan (1887-1920), native to India, 
was an extraordinary mathematician. A child prodigy, he was largely self- 
educated. When he was 16, he found an 1856 book by G. S. Carr [C], that 
listed theorems and formulas and some short proofs. Using this tutorial 
textbook, packed with facts from advanced calculus, geometry, and classical 
analysis, as a guide, Ramanujan taught himself mathematics, and by the 
age of 17 was engaged in deep mathematical research, studying Bernoulli 
numbers and divergent series and calculating the Euler-Mascheroni constant 
to 15 decimal places. 

After several unsuccessful attempts to have his work appreciated by other 
mathematicians, he wrote to G. H. Hardy, who recognized his genius and 
invited him to study with him at Cambridge. Ramanujan, a devout Brahmin, 
at first refused to travel to a foreign country, but relented when the family 
goddess Namagiri appeared to his mother in a dream commanding her not 
to prevent his departure. 

Ramanujan was tutored by, and collaborated with, Hardy for almost five 
years beginning in 1914. They published seven joint papers in varous journals 
|Ram2| . Their brilliant guest made a deep impression on Hardy and Little- 
wood, who compared him to Jacobi and Euler. Hardy considered that his 
most important mathematical achievement was the discovery of Ramanujan. 

Ramanujan was in poor health and was hospitalized for a long time be- 
cause of a non-diagnosed illness. He returned in poor health to India in 1919, 
and died soon after at the age of 32. A modern analysis of his medical records 
has indicated that he may have been suffering from a form of hepatitis. 

1.2. Ramanujan's mathematical heritage. It is beyond our competence 
to evaluate the significance of Ramanujan's mathematical genius. In the 
literature, he is sometimes mentioned in the company of other great math- 
ematicians such as Gauss, Jacobi, and Euler. Ramanujan recorded most of 
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his work in notebooks containing thousands of mathematical formulas and 
results. In a series of books, published by Springer- Verlag from 1985 through 
1998 [Belj- |Be5] , B. Berndt carefully analyzed these notebooks, giving proofs 
for the results that Ramanujan had stated without proof. Other work by 
Ramanujan is contained in the so-called Lost Notebook [Rami] and in some 
loose papers. 

Berndt's analysis and book-writing project has required an enormous 
amount of effort and scientific detective work, for which he has received 
grateful acknowledgment from the mathematical community. In particular, 
in 1996 he received the Steele Prize for Mathematical Exposition, with a ci- 
tation |S96| that reads, in part, "In an impressive scholarly accomplishment 
spread out over 20 years, Berndt has provided a readable and complete ac- 
count of the notebooks, making them accessible to other mathematicians. 
Ramanujan's enigmatic, unproved formulas are now readily available, to- 
gether with context and explication, often after the most intense and clever 
research efforts on Berndt's part." 

During the past ninety years that have passed since Ramanujan's death, 
his influence on several areas of mathematics such as number theory, com- 
binatorics, and mathematical analysis has been significant and continues to 
be so. 

1.3. AVV meets Ramanujan. For about twenty-five years the present au- 
thors had an active collaboration with the late M. K. Vamanamurthy, who 
died in 2009. We referred to our research group as AVV, after the initials 
of our last names. Our joint research dealt with geometric function the- 
ory, more precisely quasiconformal mapping theory. An important aspect of 
our work dealt with conformal invariants, usually expressed in terms of spe- 
cial functions such as the Euler gamma function, hypergeometric functions, 
complete elliptic integrals, and elliptic functions. 

By a lucky chance we discovered the survey of Askey |As] and Berndt's 
series of books on Ramanujan's notebooks [BelJ-[Be5j, where we found valu- 
able pieces of information. What interested us most was Ramanujan's work 
on the gamma and hypergeometric functions and modular equations. We 
also had access to a preprint version of [BcBGJ, which dealt with some of 
Ramanujan's theories. 

These results of Ramanujan seemed to fit nicely into our AVV research 
program, started in about 1984-in particular, the part that dealt with special 
functions, |AVV1| - |AVV5] and jABRVV] . [PVT] . |AVV6] . |AQVV| . 

1.4. Old and new research. We have written several surveys on our AVV 
research. In [AV V1| we provided an overview of the known results, along 
with several new ones, and formulated a long list of open problems, including 
problems about the gamma function. Next, in |AVV5| . we outlined some of 
our earlier results and suggested that known inequalities and identities for 
the complete elliptic integral 

%) = k(l.!;^ 2 ) 
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might hold for 2-F1 (a, b; a+b; r 2 ) with (a, b) close to (i, |) , where 2-^1(0, 6; «+ 
6; r 2 ) stands for the Gaussian hypergeometric function (see |AS| . |QLBC] and 
ff47T]) below). 

It turned out that some of these ideas bore fruit, perhaps more than 
we had anticipated, in the following years, and the results are surveyed in 
|AVV7| and [AVV8]. Several of our research topics that were inspired by 
Ramanujan's work are the Euler-Mascheroni constant, the Euler gamma 
function, volumes of unit balls in euclidean n-space, approximation of the 
Gaussian hypergeometric function 2F1, approximation of the perimeter of 
an ellipse, and the study of generalized modular equations. Some of the 
most recent results, motivated by AVV work, include the papers authored 
by X. Zhang, G. Wang, V. Heikkala, A. Baricz, E. A. Karatsuba, H. Alzer, 
and others. 

The surveys |AVV1| and |AVV5| were written before our book [AVV6J, 
which summarizes most of our work, whereas the surveys |AVV7| and |AVV8| 
were written after the publication of |AVV6| . The purpose of the present 
survey is to give a modified version of |AVV7] . with an attempt to provide 
an overview of the most recent work on this subject matter. 

While we where working on the book [AVV6], we became acquainted 
with the work of S.-L. Qiu, who subsequently visited each of the AVV team 
members at our respective home universities. He helped us in checking the 
early versions of the book manuscript, and our collaboration with him led 
to many significant co-authored results in joint papers ( |AnQ| , |AQVV| , 
|QVa| , |QVuf I , |QVu2| , |QVu3| ). Since then, his students and co-authors 
have energetically investigated the problems left open by our work, and also 



contributed in other ways to this area of research (|QF|, |QH|, |QM|, |QZ|, 
[WZG] , |WZQG| , |ZWG1IIZWG2| V 

2. Gamma function and Euler-Mascheroni constant 
Throughout this paper V will denote Euler's gamma function, defined by 



T{z) = / e~H z " l dt : Rez > 0, 
Jo 



and then continued analytically to the finite complex plane minus the set 
of nonpositive integers. The recurrence formula T(z + 1) = zT(z) yields 
r(n+ 1) = n\ for any positive integer n. We also use the fact that T(l/2) = 
y/ir. The beta function is related to the gamma function by B(a,b) = 
r(a)r(6)/r(a + b). The logarithmic derivative of the gamma function will 
be denoted, as usual, by 

The Euler-Mascheroni constant 7 is defined as (see |A1| . |TY| . |Yoj ) 

n 

7 = lim D n = 0.5772156649 . . . ; D n = V - - log n. 

n— >oo — J k 

fc=l 

The convergence of the sequence D n to 7 is very slow (the speed of conver- 
gence is studied by Alzer [All ). D. W. DeTemple |De| studied a modified 
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sequence which converges faster and proved 



1 1 n 1 / 1\ 

<jRn_7< 24^' Where Rn =^k~ log { n+ 2) 

i— 1 \ / 



24(n + l) 2 
Now let 

h(n) = R n — 7, H(n) = n 2 h(n), n ^ 1. 

Since ^(n) = —7 — 1/n + ^fc=i we see that 

H{n) = (Rn - 7 )n 2 = f *(n) + i - log (n + J n 2 . 

Some computer experiments led M. Vuorinen to conjecture that iJ(n) in- 
creases on the interval [1, 00) from H(l) = —7 + 1 — log(3/2) = 0.0173 . . . 
to 1/24 = 0.0416.... E. A. Karatsuba proved in [Kl] that for all integers 
n ^ l,H(n) < H(n + 1), by clever use of Stirling's formula and Fourier 
series. Moreover, using the relation 7 = 1 — r'(2) she obtained, for k ^ 1, 

12fe+l 12Jfc+l j f , s 

-c k ^j-l + {logk) d(k,r)- 7^Y^ c fc' 
r=X r=l 

where 

k r+1 



c k = 7Tir^ + 2k2e ~^ d(k,r) = (-l) r ' 1 



(12k)\ » v w v y ( r _i)|( r + i)' 

giving exponential convergence. Some computer experiments also seemed to 
indicate that ((n + l)/n) 2 H{n) is a decreasing convex function. 

2.1. In "The Lost Notebook and Other Unpublished Papers" of Ramanujan 
[Rami J, appears the following record: 

•T(l + ,) = ^(g'{8,3 + 4^ + , + |} I/6 , 

where X is a positive proper fraction 

30 

9 = — = .9675 

IT 6 

1/12 = .8071 7 /i2 = -3058 

02/ia = .6160 8/12 = .3014 

3/ i2 = .486 7 9/ i2 = .3041 

04/i 2 = .4029 0io/i 2 = .3118 

^5/12 = -3509 0n /12 = .3227 

6/12 = .3207 0i = .3359 

600 = 1-" 

Of course, the values in the above table, except 0^, are irrational and 
hence the decimals should be nonterminating as well as nonrecurring. The 
record stated above has been the subject of intense investigations and is 
reviewed in [BCK, page 48, Question 754]. This note of Ramanujan led the 
authors of |AVV6j to make the following conjecture. 
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2.2. Conjecture. Let 

G{x) = {e/x) x T{l + x)/V^ 

and 

H(x) = G(xf - 8x 3 - 4x 2 - x = — . 
y ' y 1 30 

Then H is increasing from (l,oo) into (1/100,1/30) |AVV6I p. 476]. 



2.3. In a nice piece of work, E. A. Karatsuba |K2| proved Conjecture 12.21 
She did so by representing the function H (x) as an integral for which she was 
able to find an asymptotic development. Her work also led to an interesting 
asymptotic formula for the gamma function: 

(2.4) 

^(x\*( , , 2 1 11 79 3539 



eJ V 30 240x 3360x 2 201600x 3 

9511 10051 47474887 07 On A 1 /® 

~ 403200x 4 ~ 716800x 5 + 1277337600x 6 + X? + " ' + tf» + n+1 ^J ' 
where A n+ i(x) = 0( n+i ), as x — )• 00, and where each is given explicitly 
in terms of the Bernoulli numbers. 



G. Nemes has studied the Ramanujan-Karatsuba formula in (|2.4p and 
shown that it is better than some other well-known approximations for the 
gamma function [N] , 

The Monotone l'Hopital's Rule, stated in the next paragraph, played an 
important role in our work [AVV4]- |AVV5] . The authors discovered this 
result in |AVV4] . unaware that it had been used earlier (without the name) 
as a technical tool in differential geometry. See |Cha| p. 124, Lemma 3.1] or 
|AQVV, p. 14] for relevant remarks. 

2.5. Lemma. For —00 < a < b < 00, let g and h be real-valued functions 
that are continuous on [a, b] and differentiate on (a, b), with U 7^ on (a, b). 
Ifg'/h' is strictly increasing {resp. decreasing) on (a, b), then the functions 

g(x) - g(a) g(x) - g(b) 

h(x) - h(a) an h(x) - h(b) 
are also strictly increasing (resp. decreasing) on (a, b). 

2.6. Monotonicity properties. In |AnQ| it is proved that the function 

(2.7) m a '^±21 

x logx 

is strictly increasing from (1, 00) onto (1 — 7, 1). In particular, for x £ (1, 00), 

(2.8) x (i-7)*-i < r(x ) < x x-i 

The proof required the following two technical lemmas, among others: 
2.9. Lemma. The function 

00 

1 \ n — x 

g{x) = > 

*— J [n + xr 

n=l K ' 

is positive for x G [1,4). 
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2.10. Lemma. The function 

(2.11) h(x) = x 2 f'(l + x) - x + x) + logT(l + x) 

is positive for all x E [1, oo). 

It was conjectured in |AnQ| that the function / in (|2,7p is concave on 
(l,oo). 

2.12. Horst Alzer |A1] has given an elegant proof of the monotonicity of 
the function / in (|2.7p by using the Monotone l'Hopital's Rule 12.51 and the 
convolution theorem for Laplace transforms. In a later paper [A3] he has 
improved the estimates in (|2.8p to 

(2.13) < r(x) < a ^(*-l)-7 ) x € (o, 1), 

where a = 1 - 7 = 0.42278 . . . , /3 = | (vr 2 /6 - 7) = 0.53385 . . . are best 
possible. If x E (l,oo), he also showed that (|2,13p holds with best constants 



a 



I (J2 



vr 2 /6 - 7) = 0.53385 (3 = 1. 



2 

2.14. Elbert and Laforgia [ELJ have shown that the function g in Lemma [279 
is positive for all x > — 1. They used this result to prove that the function 
h in Lemma 12.101 is strictly decreasing from (—1,0] onto [0, 00) and strictly 
increasing from [0,oo) onto [0, 00). They also showed that /" < for x > 1, 
thus proving the Anderson-Qiu conjecture |AnQ| , where / is as in (|2.7p . 

2.15. Berg and Pedersen [BPlJ have shown that the function / in (|2.7p is not 
only strictly increasing from (0, 00) onto (0,1), but is even a (nonconstant) 
so-called Bernstein function. That is, / > and /' is completely monotonic, 
i.e., /' > 0, /" < 0, > 0, . . . . In particular, the function / is strictly 
increasing and strictly concave on (0, 00). 

In fact, they have proved the stronger result that 1// is a Stieltjes trans- 
form, that is, can be written in the form 

1 f°° da(t) 

- c+ / — ^r, x > 0, 



f{x) J x + t' 

where the constant c is non-negative and a is a non-negative measure on 
[0, 00) satisfying 

da{t) 



< 00. 



1 + t 

In particular, for 1// they have shown by using Stirling's formula that c = 1. 
Also they have obtained da(t) = H(t)dt, where H is the continuous density 



H(t) 



iog|r(i-*)| + (fc-i)iogt t£(k _ lk)k _ 12 

(log |r(l-t)|) 2 + (fc-l) 2 7T 2 ' K — 1, 2, . 

, t = 1,2, 



Here log denotes the usual natural logarithm. The density H(t) tends to I/7 
as t tends to 0, and a has no mass at 0. 

2.16. Remark. In a series of papers I. Pinelis (see, e.g. [F]) has advocated the 
use of the Monotone l'Hopital's Rule, Lemma 12.51 Probably partly because 
of his work, during the past few years this result has found numerous appli- 
cations to the study of special functions. In a forthcoming paper |KVV] . a 
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long list of papers is provided in which Jordan's inequality is refined. Most 
of these refinements make use of Lemma 12.51 



3. Volumes of balls 

Formulas for geometric objects, such as volumes of solids and arc lengths 
of curves, often involve special functions. For example, if Q n denotes the 
volume of the unit ball B n = {x : |x| < 1} in M n , and if w n _i denotes the 
(n — l)-dimensional surface area of the unit sphere S 1 ™ -1 = {x : \x\ = 1}, 
n ^ 2, then 

Qn ~ r((n/2) + 1) ; " nfV 
It is well known that for n 7 both Vt n and co n decrease to (cf. |AVV6| 
2.28]). However, neither Q n nor uo n is monotone for n on [2, oo). On the other 
hand, Qn^ nlogn ^ decreases to e -1 / 2 as n — > oo |AVVH Lemma 2.40(2)]. 

In 2000 H. Alzer |A2| obtained the best possible constants a, b, A, B, a, 
(3 such that 

n n 

a O n+1 < O < hQ n+1 



n + A £l n -i , jn + B 

5> — 



2tt f2 n V 2tt 



1 + - < o o" ^ 1 + - 

n) iZ n _iiZ n+ i \ n 

for all integers n ^ 1. He showed that a = = 1.12837. . ., b = ^[e = 

1.64872... , A = 1/2, B = tt/2 - 1 = 0.57079... , o = 2 - (log7r)/log2 = 
0.34850 . . . , P = 1/2. For some related results, see jKiR] . 

Recently H. Alzer |A3] has obtained several sharp inequalities for f2 n . In 
particular, he showed that 

A , ..Mn+l Q n ^ B 

— 1= ^ [n + 1)— n— < — =, for n ^ 2, 



with the best possible constants A = (4 — vr)\/2 = 1.2139... and B = 
^V^tt = 1.2533 . . ., refining and complementing work in [K1RJ. 

The most recent studies dealing with the monotonicity properties of Q r 
include the following papers: |A4], [BP2] . [Ml] . [QG]. 



4. Hypergeometric functions 

Given complex numbers a, b, and c with c ^ 0, —1, —2, . . . , the Gaussian 
hypergeometric function is the analytic continuation to the slit plane C \ 
[1, oo) of 

(4.1) F(a, b; c; z) = 2 F x (a, b; q z)s ( °' " } \z\ < 1. 

jj=0 (c, n) n! 

Here (a, 0) = 1 for and (a, n) is the shifted factorial function 

(a, n) ee a(a + l)(a + 2) • • • (a + n - 1) 

for n = 1, 2, 3, . . .. 
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The hypergeometric function w = F(a, b; c; z) in (|4.ip has the simple 
differentiation formula 

(4.2) 4- F ( a > ^ c;z) = —F(a + l,b + l;c + 1; z). 

dz c 

The behavior of the hypergeometric function near z = 1 in the three cases 
a + b < c, a + b = c, and a + b > c, a, 6, c > 0, is given by 



(4.3) 



(F(a b-c-D - r (c)r(c-°-6) , fc< 

B(a, bji^a, 6; a + 6; z) + log(l - z) 

= R(a,b) + 0({l-z)log(l-z)), 
F(a, 6; c; z) = (1 - z) c ~ a ~ b F(c - a,c - b;c; z), c < a + b, 



where 

(4.4) R(a, b) = -27- *(o) - ¥(&), fl(a) = 12(a,l - 0), i2(-) = logl6, 

and where log denotes the principal branch of the complex logarithm. The 
above asymptotic formula for the zero-balanced case + 6 = c is due to 
Ramanujan (see |As| . |Bel] ). This formula is implied by [AS, 15.3.10]. 

The asymptotic formula (|4.3p gives a precise description of the behavior 
of the function F(a, b;a + b; z) near the logarithmic singularity z = 1. This 
singularity can be removed by an exponential change of variables, and the 
transformed function will be nearly linear. 



In |QF| it is shown that f(x) = R(x)sm(irx) is decreasing from (0,1/2] 
onto (tt, log 16], where the Ramanujan constant R{x) is as defined in fj4.4[> . 

4.5. Theorem. |AQVV| For a,b > 0, let k(x) = F(a,b;a + b;l-e- x ), x > 0. 
Then k is an increasing and convex function with k'((0, 00)) = (ab/(a + b), 
T{a + b)/{T{a)T{b)))- ' 

4.6. Theorem. |AQVV| Given a,b > 0, and a + b > c, d = a + b — c, the 

function £(x) = F(a,b;c; 1 — (1 + x)~ l l d ), x > 0, is increasing and convex, 
with f((0,oo)) = (ab/(cd), T (c)T (d) / (T (a)T (b))) . 

4.7. Gauss contiguous relations and derivative formula. The six func- 
tions F(a±l, b; c; z), F(a, 6±1; c; z), F(a, b; c±l; z) are said to be contiguous 
to F(a, b; c; z). Gauss discovered 15 relations between F{a, b; c; z) and pairs 
of its contiguous functions |AS| 15.2.10-15.2.27], [Rai2, Section 33]. If we 
apply these relations to the differentiation formula (|4.2p . we obtain the fol- 
lowing useful formulas. 

4.8. Theorem. For a,b,c > 0, z e (0, 1), let u = u(z) = F(a — 1, 6; c; z), 
v = v(z) = F(a, b; c; z), u\ = u(l — z), v\ = v(l — z). Then 

(4.9) z^ = (a-l)(v-u), 

az 

dv 

(4.10) z(l - z)— = (c- a)u+ (a - c + bz)v, 
and 

(4.11) —z(l - z)F(a + 1, b + 1; c + 1; z) = (c - a)u + (a - c + bz)v. 

c 
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Furthermore, 

(4.12) 2(1— z )~j£~ (uvi+mv— wi) = (1— a— b) [(1— zu\v— (1— 22)m;i )] . 

Formulas (|4,9p - (|4,lip in Theorem 14,81 are well known. See, for example, 
|AAR| 2.5.8]. On the other hand, formula (|4.12p . which follows from 
(^TU|) is first proved in |AQVV[ 3.13 (4)]. 



Note that the formula 
dF 

(4.13) 2(1 -z)— = {c- b)F(a, b - 1; c; z) + (b - c + az)F(a, 6; c; 2) 
follows from (|4,10p if we use the symmetry property F(a, b\ c; 2) = F(b, a; c; 2) 



4.14. Corollary. With the notation of Theorem KEi if a £ (0, 1), b = 1 — a < 

c, then 

(r(c)) 2 



UVl + U\V — VV\ = u(l) 



r(c + o-i)r(c-o + i) 



5. Hypergeometric differential equation 

The function F(a,b; c; 2) satisfies the hypergeometric differential equation 

(5.1) 2(1 - z)w" + [c - (a + b + l)z]w' - abw = 0. 

Kummer discovered solutions of (|5.ip in various domains, obtaining 24 in 
all; for a complete list of his solutions see |Rai2| pp. 174, 175]. 

5.2. Lemma. (1) If2c = a + b+l then both F(a, b; c; 2) and F(a, 6; c; 1 — 2) 
satisfy (|5.ip in the lens-shaped region {z : < \z\ < 1, < |1 — z\ < 1}. 
(2) If 2c = a + b + 1 i/ien 6o£/i -F(a, b\ c; 2 2 ) and .F(a, 6; c; 1 — 2 2 ) satisfy the 
differential equation 

(5.3) 2(1 - 2 2 )w" + [2c - 1 - (2a + 26 + l)2 2 ]u/ - 4abzw = 

in t/ie common part of the disk {2 : |z| < 1} and the lemniscate {2 : |1 — 2 2 | < 
1}. 

Proof. By Kummer (cf. |Rai2| pp. 174-177]), the functions F(a,b;c; z) 
and F(a, b; a + b + 1 — c; 1 — 2) are solutions of (|5.ip in {2 : < \z\ < 1} and 
{2 : < |1 — z\ < 1}, respectively. But a + 6+1 — c = c under the stated 
hypotheses. The result (2) follows from result (1) by the chain rule. □ 

5.4. Lemma. The function F(a, b; c; y/l — z 2 ) satisfies the differential equa- 
tion 

Z 3 (l - Z)zw" - {Z(l - Z) + [c-{a + b+ l)Z]Zz 2 }w' - abz 3 w = 0, 

in the subregion of the right half-plane bounded by the lemniscate r 2 = 
2cos(2tf), -71-/4 ^ ■& ^ tt/4, 2 = re il9 . Here Z = y/l - z 2 , where the 
square root indicates the principal branch. 
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Proof. From (4.1), the differential equation for w = F(a,b;c;t) is given 



by 



. d 2 w r , , . . dw 

*(! - t)-TT + [c-(a + b+ l)t]— abw = 0. 

dt z dt 



Now put t = y/l — z 2 . Then 

dz _ t dt z d 2 t 1 

dt z' dz t' dz 2 t 3 

and 

dw t dw d 2 w t 2 d 2 w 1 dw 

dt z dz' dt 2 z 2 dz 2 z 3 dz 

So 



t 2 „ 1 ,1 r , . ^ .1 / t 



i(l-t) ^w" + c- (a + b + l)t 

-Z i z 6 



— )w' — abw = 0. 

Z' 



Multiplying through by z 3 and replacing t by Z = y/l — z 2 gives the result. □ 

If w\ and W2 are two solutions of a second order differential equation, then 
their Wronskian is defined to be W(w\, w%) = wiw' 2 — W2w[. 

5.5. Lemma. |AAR| Lemma 3.2.6] Ifw\ andw2 are two linearly independent 
solutions of (|5.ip . then 

W(z) = W{ Wl ,w 2 ){z) A 



Z c (l - z )a.+b-c+l ' 

where A is a constant. 

(Note the misprint in |AAR| (3.10)], where the coefficient x(l — x) is 
missing from the first term.) 



5.6. Lemma. If 2c = a + b + 1 then, in the notation of Theorem \4.8\. 
(5.7) (c - a)(uvi + uiv) + (a - l)vvi = A ■ z x ~ c {\ - z) 1 ' . 

For a proof see |AVV7| . Note that in the particular case c = 1, a = b = ^ 
the right side of (|5.7p is constant and the result is similar to Corollary 14.141 
This particular case is Legendre's Relation (|6.3| an elegant proof of it was 
given by Duren |Du| . 

5.8. Lemma. If a, b > 0, c ^ 1, and 2c = a + b + 1, then the constant A 
in Lemma 15.61 is given by A = (T(c)) 2 /(T(a)T(b)). In particular, if c = 1 
then Lemma 15.61 reduces to Legendre 's Relation (|6.6p for generalized elliptic 
integrals. 

For a detailed proof of this lemma we refer the reader to [A VV7| . 



For rational triples (a, b, c) there are numerous cases where the hyper- 
geometric function F(a,b;c; z) reduces to a simpler function (see jPBMj). 
Other important particular cases are generalized elliptic integrals, which we 
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will now discuss. For a, r £ (0, 1), the generalized elliptic integral of the first 
kind is given by 

X a = %a(r) = ^F(a,l-a;l;r 2 ) 

rn/2 

= (sin7ra) / (tant) 1-2o (l - r 2 sin 2 t)~ a alt, 
Jo 

%' a = %' a {r)=% a {r'). 
We also define 

Mar )= s , r'= Vl-r 2 . 

2sm(7ra) JC a (r) 



For fi a ( r ) some functional inequalities are obtained in |QH| , some interest- 
ing monotonicity properties in |WZQC| , and several sharp inequalities in 
|ZWC2| . 

The invariant of the linear differential equation 
(5.9) w" + pw + qw = 0, 

where p and g are functions of z, is defined to be 

i = q-^p --£p 

(cf. [Rai2,p.9]). If w\ and ^2 are two linearly independent solutions of (|5.9p . 
then their quotient w = Wijw\ satisfies the differential equation 

S w (z) = 21, 

where S w is the Schwarzian derivative 



and the primes indicate differentiations (cf. [Rai2, pp. 18,19]). 

From these considerations and the fact that % a { r ) an d 9C' a (r) are linearly 
independent solutions of (|5.3p (see |AQVV (1.11)]), it follows that w = fJ- a ( r ) 
satisfies the differential equation 

-8a(l - a) 1 + 6r 2 - 3r 4 
<WJ - + 2r 2 (r') 4 ' 

The generalized elliptic integral of the second kind is given by 

£a = £a(r) = ^F(a- 1,1- a; l;r 2 ) 
/■tt/2 

= (sin™) / (tan t) 1 " 20 ^ - r 2 sin 2 t) 1_a dt 
£' a = £' a (r) = £«(/), 

£.(o)=;, £.(i) - sin(7ra) 



2' " w 2(1 - a)' 

For a = |, 3C a and £ a reduce to 3C and £, respectively, the usual elliptic 
integrals of the first and second kind [BFJ, respectively. Likewise [X\ji (r) = 
n{r), the modulus of the well-known Grotzsch ring in the plane [LVJ. 
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5.10. Corollary. The generalized elliptic integrals % a and £ a satisfy the dif- 
ferential equations 

(5.11) r ^ r ') 2 ^r + C 1 - 3r " 2 )^ - 4a (! " a > % - = °> 

(5.12) K^') 2 ^ + (r'f^- + 4(1 - afrE a = 0, 
respeciiwe/y. 

Proof. These follow from JO}. □ 

For a = o these reduce to well-known differential equations |AVV6( pp. 
474-475], [BF] . 



6. Identities of Legendre and Elliott 

In geometric function theory the complete elliptic integrals %{r) and £(r) 
play an important role. These integrals may be defined, respectively, as 

%{r) = f F(i, I; 1; r 2 ), £(r) = f -±; 1; r 2 ), 

for — 1 < r < 1. These are 3C a (r) and £ a (?~), respectively, with a = ^. We 
also consider the functions 

%' = %'{r) = %( r '), < r < 1, 

3C(0) = 7r/2, 3C(l - ) = +oo, 

and 

£' = £'(r) = £(r'), O^r^l, 

where r' = \/l — f 2 - For example, these functions occur in the following 
quasiconformal counterpart of the Schwarz Lemma |LV] : 



6.1. Theorem. For K G [l,oo), let w be a K -quasiconformal mapping of 
the unit disk D = {z : \z\ < 1} into the unit disk D' = {w : \w\ < 1} with 
w (0) = 0. Then 

\w(z)\ ^ if K (\z\), 

where 

(6.2) ^( r ) = M -i(l M ( r )) an d M ( r ) = |£M. 

77ms result is sharp in the sense that for each z G -D anc? G [1, oo) i/iere is 
an extremal K -quasiconformal mapping that takes the unit disk D onto the 
unit disk D' with u>(0) = and \w(z)\ = (px(\z\) (see |LV| p. 63]). 

It is well known [BF] that the complete elliptic integrals X and £ satisfy 
the Legendre relation 

(6.3) EX' + E'X-XX' = |. 

For several proofs of f j 6 . 3 [) see |Du] , 

In 1904, E. B. Elliott [E] (cf. |AVV5] ) obtained the following generalization 
of this result. 



Reflections on Ramanujan's Mathematical Gems 



13 



6.4. Theorem. If a,b, c ^ and < x < 1 then 

ffi ^ /? p _i_ /? p pp_ r(g + 6+i)r(6 + c +i) 

r(a + 6 + c+ |)r(6+ ~, 

where 

Fl = i? (^ + a '~^~ C;1 + a + 6 ' X )' 
i ? 2 = i ? Q-a,^ + c;l + 6 + c;l 

i ? 3 = i ? Q + a,^-c;l + a + 6;x 

F 4 = i*M - - - a, - + c; 1 + b + c; 1 - .>■ 



Clearly (|6.3p is a special case of (|6.5p . when a = 5 = c = and x = 
r 2 . For a discussion of generalizations of Legendre's Relation see Karatsuba 
and Vuorinen |KV] and Balasubramanian, Ponnusamy, Sunanda Naik, and 
Vuorinen |BPSV] , 

Elliott proved fj6. 5 [) by a clever change of variables in multiple integrals. 
Another proof was suggested without details in |AAR| p. 138], and in |AVV7| 
we provided the missing details. 

The generalized elliptic integrals satisfy the identity 

(6-6) S. a %' a + E' a X a - ^a^'a = Ar ', - — r- 

4(1 - a) 

This follows from Elliott's formula (|6.5p and contains the classical relation 
of Legendre (|6.3p as a special case. See also Lemma 15.81 

Finally, we record the following formula of Kummer [Kuml p. 63, Form. 
30]: 

F(a, b;a + b - c + l;l - x)F(a + 1, b + 1; c + 1; x) 

-F(a, b; c; x)F(a + l,6+l;a + 6- c + 2;l-x) 



a+b-c+1 

--Dx- c {l-xf- a -°-\ D 



c „ a _ b _x ^_T(a + b-c + l)T(c+l) 



r(o + i)r(6 + i) 

This formula, like Elliott's identity, may be rewritten in many different ways 
if we use the contiguous relations of Gauss. Note also the special case c = 
a + b-c+1. 



7. Approximation of elliptic integrals and perimeter of 

ELLIPSE 

Efficient algorithms for the numerical evaluation of %{r) and £(r) are 
based on the arithmetic-geometric mean iteration of Gauss. This fact led to 
some close major ant /minor ant functions for 0C(r) in terms of mean values 
in |VV| . Recently, mean iterations derived from transformation formulas for 
the hypergeometric functons have been investigated in jHKMj. 
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Next, let a and b be the semiaxes of an ellipse with a > b and eccentricity 
e = \J o? — b 2 /a, and let L(a, b) denote the perimeter of the ellipse. Without 
loss of generality we take a = 1. In 1742, Maclaurin (cf. [ABJ) determined 
that 

L(l,6) = 46(e) = 27r- 2 Fi(i, -±;l;e 2 ). 
In 1883, Muir (cf. |AB| ) proposed that L(l,6) could be approximated by 
the expression 27r[(l + 6 3 / 2 )/2] 2 / 3 . Since this expression has a close resem- 
blance to the power mean values studied in |VV] . it is natural to study the 
sharpness of this approximation. Close numerical examination of the error 
in this approximation led Vuorinen [V] to conjecture that Muir's approxi- 
mation is a lower bound for the perimeter. Letting r = \f\ — 6 2 , Vuorinen 
asked whether 

C7.D ^^(i-i^^^r 

for all r E [0,1]. 

In [BPRlJ Barnard and his coauthors proved that inequality (|7.ip is true. 
In fact, they expanded both functions into Maclaurin series and proved that 
the differences of the corresponding coefficients of the two series all have the 
same sign. 

Later, the same authors [BPR2J discovered an upper bound for £ that 
complements the lower bound in (j7.ip : 

(7-2) ^£(r) = 2 Fi(i-i;l;r 2 )<(i±|^) 1/2 , (Ur<l. 

See a lso (BPS ]. 

In [BPR2J the authors have considered 13 historical approximations (by 
Kepler, Euler, Peano, Muir, Ramanujan, and others) for the perimeter of an 
ellipse and determined a linear ordering among them. Their main tool was 
the following Lemma 17.31 on generalized hypergeometric functions. These 
functions are defined by the formula 

i(aj,n) z 11 

p F q (ai,a 2 ,--- ,a p ;bi,b 2 ,--- ,b q ;z) = 1 + q • — , 

B =l U j=lH' n J U - 

where p and q are positive integers and in which no denominator parameter 
bj is permitted to be zero or a negative integer. When p = 2 and q = 1, this 
reduces to the usual Gaussian hypergeometric function F(a, b; c; z). 
Some of this joint research is discussed in the survey paper |BRT| . 

7.3. Lemma. Suppose a,b > 0. Then for any e satisfying oj^o; < e < 1, 

3 F 2 (-n, a, b; 1 + a + b, 1 + e - n; 1) > 
/or a// integers n ^ 1. 

7.4. Some inequalities for 3C(r). At the end of the preceding section we 
pointed out that upper and lower bounds can be found for %(r) in terms 
of mean values. Another source for the approximation of %(r) is based on 
the asymptotic behavior at the singularity r = 1, where %{r) has logarith- 
mic growth. Some of the approximations motivated by this aspect will be 
discussed next. 
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Anderson, Vamanamurthy, and Vuorinen |AVV3| approximated 3C(r) by 
the inverse hyperbolic tangent function arth, obtaining the inequalities 

/ \ 1 I 2 

,„ 7r / arthr \ „, . n arthr 

< 75 » 2 {—) <*(')< J—- 

for < r < 1. Further results were proved by Laforgia and Sismondi [L 
Kiihnau [Kulj and Qiu [Q] proved that, for < r < 1, 

9 < X(r) 



8 + r 2 log(4/r')' 

Qiu and Vamanamurthy |QVa| proved that 

%(r) 1 

W < 1 + -(r'f for < r < 1. 



log(4/r') 4 

Several inequalities for %{r) are given in |AVV6| Theorem 3.21]. Later Alzer 
3] showed that 



V 4 log 2 J { ' log(4/r')' 

for < r < 1. He also showed that the constants 4 and 7r/(41og2) — 1 in 
the above inequalities are best possible. The authoritative NIST handbook 
[OLBCJ lists some of these inequalities in its Section 19.9. 
For further refinements, see |QVul[ (2.24)] and |Be| . 

Alzer and Qiu |A1Q| have written a related manuscript in which, besides 
proving many inequalities for complete elliptic integrals, they have refined 
(|7.5p by proving that 

7T/arthr\ 3 / 4 , , ir arthr 
-(— ) <X(r)<- — . 

They also showed that 3/4 and 1 are the best exponents for (arth r)/r on the 
left and right, respectively. Further estimates for complete elliptic integrals 
have been obtained in [ABa| and |GQ| . 

One of the interesting tools used by the authors of |A1Q| is the following 
lemma of Biernaki and Krzyz [BKJ (for a detailed proof see |PV1] ): 

7.6. Lemma. Let r n and s n , n = 1, 2, . . . be real numbers, and let the power 
series R(x) = Y^=i r n xn an d S(x) = Y^=\ s nX n be convergent for \x\ < 
1. If s n > for n = 1,2,..., and if r n /s n is strictly increasing {resp. 
decreasing) for n = 1,2,..., then the function R/S is strictly increasing 
{resp. decreasing) on (0,1). 

7.7. Generalized elliptic integrals. Recently some new estimates for %{r) 
and £(r) were obtained in |GQ| . For the case of generalized elliptic integrals 
some inequalities are given in |AQVV , and further properties are found for 



them in [HLVV] and |HVV| . B. C. Carlson has introduced some standard 
forms for elliptic integrals involving certain symmetric integrals. Approxima- 
tions for these functions can be found in [CGJ. In |KN] and |OLBC| several 
inequalities are obtained for elliptic integrals given in the Carlson form. In 
[ABaJ S. Andras and A. Baricz have compared the generalized elliptic in- 
tegral % a ( r ) with certain other zero-balanced hypergeometric functions. In 



16 



G. D. ANDERSON AND M. VUORINEN 



his new book Baricz |Bar3| investigates various properties of power series 
and provides refinements for some of the above results, applying, for ex- 
ample, Lemma 17.61 See also Zhang, Wang and Chu [ZWClJ. Recently, a 
variant of Lemma 17.61 for the case when the numerator and denominator are 
polynomials of the same degree, was give in |HVV] . See also [KSJ. 

8. Hypergeometric series as an analytic function 

For rational triples (a, b, c) the hypergeometric function often can be ex- 
pressed in terms of elementary functions. Long lists with such triples con- 
taining hundreds of functions can be found in |PBM| . For example, the 
functions 

f(z) = zF(l 1 l;2;z) = -log(l-z) 

and 

g(z) = zF(l,l/2;3/2;z 2 ) = ^log 

have the property that they both map the unit disk into a strip domain. 
Observing that they both correspond to the case c = a + b one may ask 
(see |PV1| IPV2] ) whether there exists 5 > such that zF(a, b;a + b; z) and 
zF(a, b;a + b; z 2 ) with a, b E (0, 5) map into a strip domain. 

Membership of hypergeometric functions in some special classes of univa- 
lent functions is studied in |PVT1 lPV2l lPV3l IBPV2] . 

9. Generalized modular equations 

The argument r is sometimes called the modulus of the elliptic integral 
%(r); further, for integer values p = 1, 2, . . . , the equation 

(9.D %&=* {r) 



%(s) * X(r) ' 

with r, s E (0,1), is called the modular equation of degree p. If we use the 
notation 

(9.2) VK {r) = ^\^{K)/K), 

where \x is the modulus of the well-known Grotzsch ring in the plane 11. Vj . 
then the solution of (|9.ip for s is given by s = (fii/ p (r). We will now discuss 
some of the numerous modular equations or, more precisely, algebraic conse- 
quences of the transcendental equation (|9.ip that were found by Ramanujan. 
Our discussion is based on a nice survey of Ramanujan's work in |Be3| pp. 
4-10]. 

Ramanujan introduced the convenient notation 

a = r 2 , p = s 2 

for use in connection with (|9.ip . In this notation a third-degree modular 
equation due to Legendre |BB| p. 105] takes the form 

M )l/4 + ((1 _ a)(1 _ /3)) l/4 = 1) 

with a = r 2 , (3 = ^i/^{r) 2 . In the next theorem we list some of Ramanujan's 
modular equations, in his notation. 

9.3. Theorem. The function (fx satisfies the following identities for r E 
(0,1): 
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(1) {a?) 1 ' 2 + (((i _ a )(i _ 0))i/2 + 2 (16a/3(l - a)(l - /3))V6 = i 
for a = r 2 ,/3 = (pi/ 5 (r) 2 . 

(2) (a/3) 1 /8 + (( i_ a)( i_ /3)) i/8 = 1 

for a = r 2 , f3 = (fip{r) 2 . 

(3) (a(l - 7 )) 1/8 + (7(1 - a)) 1/8 = 2 1 /3( /3 (i _ /3))i/24 
/or a = r 2 ,/3 = ^ 1/3 (r) 2 ,7 = ^i/ 9 (r) 2 . 

(4) (a/3) 1 / 8 + ((1 - a)(l - /3)) 1 / 8 + 2 2 / 3 (a/3(l - a)(l - /3)) 1 / 24 = 1 
for a = r 2 , /3 = (p 1/2 3{r) 2 - 

(5) (I(l + ^+ v /(l-«)(l-/3))) 1 / 2 = ( a/3 )i/8 + ((l- a )(l-/3))i/8- 
(a/3(l-a)(l-/3)) 1 /8 

/or a = r 2 ,/3 = </?i/ 7 (r) 2 or for a = (/? 1/3 (r) 2 , /3 = ^i/ 5 (r) 2 . 

All of these identities are from |Be3j: (1) is Entry 13 (i) on p. 280, (2) is 
Entry 19 (i) on p. 314, (3) is Entry 3 (vi) on p. 352, (4) is Entry 15 (i) on 
p. 411, and (5) is Entry 21 (i) on p. 435. 

In 1995 B. Berndt, S. Bhargava, and F. Garvan published an important 
paper |BeBG| in which they studied generalized modular equations and 
gave proofs for numerous statements concerning these equations made by 
Ramanujan in his unpublished notebooks. No record of Ramanujan's origi- 
nal proofs has remained. A generalized modular equation with signature 1/a 
and order (or degree) p is 

(94) *-(q,i-q;i;i-* 2 ) _ /M-«; i;i-r 2 ) 

1 ' F(a,l-a;l;s 2 ) P F(a, 1 - a; 1; r 2 ) ' < < 



Here F is the Gaussian hypergeometric function defined in (|4.ip . The word 
generalized alludes to the fact that the parameter a S (0, 1) is arbitrary. 
In the classical case, a = \ and p is a positive integer. Modular equations 
were studied extensively by Ramanujan, see |BeBG| . who also gave numerous 
algebraic identities for the solutions s of (|9.4p for some rational values of a 

Slicll QbS g , ^ , . 

To rewrite (|9.4p in a slightly shorter form, we use the decreasing homeo- 
morphism /j, a : (0, 1) — > (0, oo) defined by 

7r F(a, 1 — a; 1; 1 — r 2 ) 
MaW = 2sin(vra) F(a, 1 - a; 1; r 2 ) ' 
for a 6 (0, 1). We can now rewrite (|9.4p as 

(9.5) fi a (s) =PHa(r), < r < 1. 
The solution of (|9.5p is then given by 

(9.6) s = <p a K (r) = »- l (» a (r)/K), p = l/K. 

We call (p%;(r) the modular function with signature 1/a and degree p = l/K. 
Monotonicity and convexity properties of % a (r), £ a (0, c / 3 A'( r )j an d Hair) and 
certain combinations of these special functions are established in |WZC| . 

For the parameter K = 1/p with p a small positive integer, the function 
(|9.6p satisfies several algebraic identities. The main cases studied in |BeBG| 
are 

1 1 1 

a = 7'7'o' P = 2,3,5,7,11,. .. . 
b 4 3 



18 



G. D. ANDERSON AND M. VUORINEN 



For generalized modular equations we use the Ramanujan notation: 



a^r 2 , P = vl /P {r?. 

We next state a few of the numerous identities |BeBG| satisfied by <4>\j v 
for various values of the parameters a and p. 

9.7. Theorem. ([BeBUJ Theorem 7.1]) // (3 has degree 2 in the theory of 
signature 3, then, with a = |, a = r 2 , (3 = ipy 2 (r) 2 , 

(a/3) l/2 + {(1 _ a)(1 _ /3)} l/3 = L 

9.8. Theorem. QBeBG] Theorem 7.6]) If (3 has degree 5 then, with a = |, 
a = r 2 , = <^ /3 (r) 2 , 

( a/3 )i/3 + {(1 _ a)(1 _ + 3{a/ g (1 _ a)(1 _ /3 )}i/6 = i 

9.9. Theorem. QBeBG] Theorem 7.8]) // /3 has degree 11 i/ien, with a = |, 
a = r 2 , /? = <^ /n (r) 2 , 

(a/3) 1 / 3 + {(1 - a)(l - /3)} 1 / 3 + 6{a/3(l - a)(l - /3)} 1/6 

+3v / 3{a/3(l-a)(l-/3)} 1 / 12 {(a/3) 1 / 6 +((l-a)(l-/3)) 1 / 6 } = 1. 

Such results are surprising, because they provide algebraic identities for 
the modular function tp a K , which itself is defined in terms of the transcen- 
dental function Hair). It is an interesting open problem to determine which 
of the modular equations in |BeBG| can be solved algebraically, explicitly in 
terms of the modular function. 

Because of its geometric significance for geometric function theory (see 
(|6.2p . |LVj . |QVu3| , |Ku2| ). it is desirable to give upper bounds for the func- 
tion ipk{t\K > 1 • There are many such bounds in the literature; see the 
survey |AVV8| . Much less is known about the function (p a K {r),K > 1 , but 
some bounds can be found in |HWj . |HLVV| . |WZC] . |WZQC| . 
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